hep-ph/0001005 
CERN-TH/2K-001 
UMN-TH-1834/2K 
TPI-MINN-2K/01 

Re-Evaluation of the Elastic Scattering of Supersymmetric Dark 

Matter 

o : 
o . 

CN ■ John Ellis\ Andrew FerstP and Keith A. Ohve^'^ 

c : 

I ^ I ^TH Division, CERN, Geneva, Switzerland 

^ ■ "^School of Physics and Astronomy, University of Minnesota, Minneapolis, MN 55455, USA 

^ ! ^Theoretical Physics Institute, University of Minnesota, Minneapolis, MN 55455, USA 

in : 

Q ■ Abstract 

o : 

O ■ We examine the cross sections for the elastic scattering of neutrahnos x nucleons p, n, 

^ ' as functions of in the constrained minimal supersymmetric standard model. We find 

P^l narrow bands of possible values of the cross section, that are considerably lower than some 

Qh! previous estimates. The constrained model is based on the minimal supergravity-inspired 

^ . framework for the MSSM, with universal scalar and gaugino masses mo, mi/2, and jj and the 

* * ■ 

■ MSSM Higgs masses treated as dependent parameters. We explore systematically the region 
^ I of the plane where LEP and other accelerator constraints are respected, and the 

I relic neutralino density lies in the range 0.1 < ^x^'^ — ^-^ preferred by cosmology. We 
update previous discussions of both the spin-independent and -dependent scattering matrix 
elements on protons and neutrons, using recent analyses of low-energy hadron experiments. 



CERN-TH/2K-001 
January 2000 



1 Introduction 



One of the most promising candidates for the cold dark matter beheved to pervade the Uni- 
verse is the hghtest supersymmetric particle (LSP) [IH, commonly expected to be the lightest 
neutralino x-, which is stable in the minimal supersymmetric extension of the Standard Model 
(MSSM) with conserved R parity 0]. The quantum stability of the gauge hierarchy sug- 
gests that sparticles weigh less than about 1 TeV 0, which is also the range favoured for a 
cold dark matter particle |Q, and there are indeed generic domains of the MSSM parameter 
space in which the relic LSP density falls within the range 0.1 < Vt^h^ < 0.3 favoured by 
astrophysics and cosmology The unsuccessful laboratory searches for sparticles impose 
non-trivial constraints on the MSSM parameter space, suggesting that the LSP x is mainly 
a f/(l) gaugino (Bino) 

Many non-accelerator strategies to search for cosmological relic neutralinos have been 
proposed including indirect searches for products of their annihilations in free space 
or inside astrophysical bodies, and direct searches for their scattering on target nuclei in 
low-background underground laboratories . The rates for such experiments typically have 
larger uncertainties than those for producing sparticles at accelerators, since they involve 
some astrophysical and/or cosmological uncertainties as well as those due to simulations of 
the signatures, over and above the common uncertainties in the MSSM parameters. Nev- 
ertheless, such dark matter searches offer interesting prospects for beating accelerators to 
the discovery of supersymmetry, particularly during the coming years before the LHC enters 
operation. 

In this paper, we embark on a programme to clarify the extents of the uncertainties 
in searches for supersymmetric dark matter, by re-evaluating the rates to be expected for 
the elastic scattering of relic LSPs on protons and neutrons Large ranges for these 

rates are often quoted reflecting general explorations of the MSSM parameter space. 
Our approach is to establish as accurately as possible a baseline set of predictions based on 
the most plausible assumptions commonly used in constrained MSSM phenomenology, such 
as universality in the soft supersymmetry-breaking parameters as suggested by minimal 
supergravity models, and requiring the cosmological relic density to lie within the range 
favoured by astrophysics and cosmology, namely 0.1 < ^x^"^ — These assumptions 

can and should be questioned, but they are well motivated and good candidates for default 
options in analyses of the MSSM and cold dark matter. 

In the course of this re-evaluation of elastic x — p, n scattering cross sections, we re- 
analyze the relevant spin-independent and spin-dependent matrix elements of scalar densities 



1 



and axial currents in protons and neutrons. We update previous analyses using further 



information from chiral symmetry [|Tl|, |T2[, low-energy n — p,n scattering |T3| and deep- 
inelastic lepton-nucleon scattering We include a discussion of uncertainties in the values 
of the scalar and axial-current matrix elements. 

We perform a systematic scan of the region of the mo, mi/2 parameter space of the 
MSSM with supergravity-inspired universality that is consistent with accelerator constraints 



and yields a cosmological relic density within the favoured range 0.1 < < 0.3 [15|. We 
treat /i as a dependent parameter (modulo a sign ambiguity), and our results are not very 
sensitive to A. We order our results in terms of rriy. which closely tracks mi/2- For any 
given choice of m^, tan /5 and the sign of /x, we find a relatively narrow band of possible cross 
sections, reflecting the fact that the accelerator and cosmological constraints 0, |15| favour 
a predominant U{1) gaugino (Bino) composition for the LSP. Our results fall considerably 



below many of the possible predictions in the literature |10|, and may discourage some faint- 
hearted experimentalists. However, we think they provide a realistic estimate of the target 
sensitivity required for an experiment to have a good chance of success. 

2 Theoretical Framework 

We review in this Section the theoretical framework we use in the context of the MSSM [g. 
The neutralino LSP is the lowest-mass eigenstate combination of the Bino B, Wino W 
and Higgsinos H12, whose mass matrix is diagonalized by a matrix Z: diag{m^-^^^^^4) = 
Z*NZ^^. The composition of the lightest neutralino may be written as 

X = Z^iB + Z^2W + Z^^Hi + Z^,H2 (1) 

As already mentioned, we assume universality at the supersymmetric GUT scale for the 
f/(l) and SU (2) gaugino masses: Mi^2 = ^^1/25 so that Mi = | tan^ 9wM2 at the electroweak 
scale. We denote by tan /? the ratio of Higgs vacuum expectation values, and yU is the Higgsino 
mass-mixing parameter. We also assume GUT-scale universality for the soft supersymmetry- 
breaking scalar masses mg, for the Higgs bosons as well as the squarks and sleptons. We 
further assume GUT-scale universality for the soft supersymmetry-breaking trilinear terms 
A. Our treatment of the sfermion mass matrices M follows [|16|]. As discussed there, the 



sfermion mass-squared matrix is diagonalized by a matrix rj: diag{mi, m^) = rjM'^r] ^, which 
can be parameterized for each flavour / by an angle 6f and phase 7/: 

(2) 
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As a simplification, we neglect CP violation in this paper, so that 7/ = and there are no 
CP-violating phases in the neutralino mass matrix, either. We treat mi/2,iTLo, A and tan (3 
as free parameters, and fi and the pseudoscalar Higgs mass rriA as dependent parameters 
specified by the electroweak vacuum conditions, which we calculate using nit = 175 GeV Q 
The MSSM Lagrangian leads to the following low-energy effective four-fermi Lagrangian 



suitable for describing elastic x-nucleon scattering [|16 



£ = xYl^XQilf^ioiii + Oi2i1^)qi + OisiXXQiQi + oiAiXl^XQil^Qi + Oi^iXXQil^ Qi + OifiiXl^Xmi (3) 

This Lagrangian is to be summed over the quark generations, and the subscript i labels up- 
type quarks {i = 1) and down-type quarks {i = 2). The terms with coefficients aii,a4i,a;5j 
and a^i make contributions to the elastic scattering cross section that are velocity-dependent, 
and may be neglected for our purposes. In fact, if the CP violating phases are absent as 



assumed here, = = |17]. The coefficients relevant for our discussion are: 



and 
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where 
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^We have checked that varying nit by ±5 GeV has a neghgible effect on our results. 
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where Ui^T^i denote hypercharge and isospin, and 

-Bj = sin/?(cos/5) , = cos/?(— sin/5), 

Cj = sina;(cosa) , = cosa(— sina) (7) 

for up (down) type quarks. We denote by rriH^ < uih^ the two scalar Higgs masses, and a 
denotes the Higgs mixing angle 0. 

3 Hadronic Matrix Elements 

The elastic cross section for scattering off a nucleus can be decomposed into a scalar (spin- 
independent) part obtained from the a2i term in (^), and a spin-dependent part obtained 
from the a^i term. Each of these can be written in terms of the cross sections for elastic 
scattering for scattering off individual nucleons, as we now review and re-evaluate. 
The scalar part of the cross section can be written as 

4?T7 ^ 

a, = ^[ZU + {A-Z)Uf (8) 

TT 

where rrir is the reduced LSP mass, 

^= E f^^i^ + y^^iY.— (9) 

P q=u,d,s y c,b,t y 

and /„ has a similar expression. The parameters f^^ are defined by 

^vfrl = {Pl^qQQlP) = ^qBq (10) 

whilst f^Q = 1 — J2q=u,d,s frq [mi- observe that only the products nigBq, the ratios of 
the quark masses rriq and the ratios of the scalar matrix elements Bg are invariant under 
renormalization and hence physical quantities. 



We take the ratios of the quark masses from |11 



in Tfi 

— = 0.553 ± 0.043, — = 18.9 ±0.8 (11) 

vfid rud 

In order to determine the ratios of the Bg and the products rrigBq we use information from 
chiral symmetry applied to baryons. Following |jl2|, we have: 



Bu — Bs m=o + m=- — m„ — m„ , , 

z=— = " — = ^ 12 

Bd - Bs mE+ + ms- - nip - rrin 



^Wc note that (||) is taken from [|T^ and corrects an error in and that (^, ||) agree with |j and 
the pubUshed version of Mm. 
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Substituting the experimental values of these baryon masses, we find 

z = 1.49 (13) 
with an experimental error that is negligible compared with others discussed below. Defining 

y ^ 75— V' 



we then have 

Bd 2 + {{z-l)xy) 



(15) 



a = i(m„ + mrf) X (5, + B„) =45±8 MeV (16) 



2xz-{{z-l)xy) 
The experimental value of the vr-nucleon a term is ||13||: 



and octet baryon mass differences may be used to estimate that []T3 



a = : ao = 36 ± 7 MeV (17) 

(1 -y) 



Comparing (|TBp and ([T^), we find a central value oi y = 0.2, to which we assign an error 
±0.1, which yields 

-^ = 0.73 ±0.02 (18) 
B„ 



The formal error in y derived from ([iq) and ([T7| ) is actually ±0.2, which would double the 
error in B^/Bu- We have chosen the smaller uncertainty because we consider a value of y in 
excess of 30% rather unlikely. However, we do illustrate later by one example the potential 
consequences of a larger error in y. 

The numerical magnitudes of the individual renormalization-invariant products rrigBq 
and hence the f^^ may now be determined: 

/jfj = 0.020 ± 0.004, /jfj = 0.026 ± 0.005, f^^] = 0.118 ± 0.062 (19) 

where essentially all the error in f^?] arises from the uncertainty in y. The corresponding 
values for the neutron are 

/^"^ = 0.014 ± 0.003, 4'^^ = 0.036 ± 0.008, /^"^ = 0.118 ± 0.062. (20) 

It is clear already that the difference between the scalar parts of the cross sections for 
scattering off protons and neutrons must be rather small. 
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The spin- dependent part of the elastic x-nucleus cross section can be written as 

32 



(72 = — G|m^AV( J + 1) 

TT 



where is again the reduced neutrahno mass, J is the spin of the nucleus, and 

1 



A = j{ap{Sp) + an{Sn)) 



where 



Q^2» ^(n) 



(21) 



(22) 



(23) 



The factors Aj-^'""* parametrize the quark spin content of the nucleon. A recent global analysis 
of QCD sum rules for the gi structure functions [Q, including 0{a^) corrections, corresponds 
formally to the values 



A^f) = 0.78 ± 0.02, aJ'^ = -0.48 ± 0.02, A'f^ = -0.15 ± 0.02 



(24) 



whilst perturbative QCD fits to the data for gi tend to give broader ranges |T^. In our 
numerical analysis, we double the formal errors in (^) to ±0.04, essentially 100% correlated 
for the three quark fiavours. In the case of the neutron, we have A^;^) = aJ^ A^"^ = A^p\ 
and Al") = Al^). 



4 Cosmological and Experimental Constraints 



The domain of MSSM parameter space that we explore in this paper is that defined in [15 



Several convergent measures of cosmological parameters ||20[ suggest that the cold dark 
matter density ^Icdm = 0.3±0.1 and that the Hubble expansion rate H = /ix 100 km/s/Mpc: 
h = 0.7±0.1, leading to the preferred range 0.1 < VtcDuh^ < 0.3. The upper limit on VLqdm 
can be translated directly into the corresponding upper limit on Vt^. However, it is possible 
that there is more than one component in the cold dark matter, so that Vl^ < Qcdm, opening 
up the possibility that < 0.1. Although the MSSM parameters which lead to < 0.1 
tend to give larger elastic scattering cross sections, the detection rate also must be reduced 
because of the corresponding reduction in the density of LSPs in the Galactic halo. Here 
we shall neglect this possibility, assuming instead that essentially all the cold dark matter is 
composed of LSPs, so that fl^ > 0.1. 



For the calculation of the relic LSP density, we follow |]I5|, where coannihilations between 
X and the sleptons i, particularly the lighter stau fi, were shown to play an important role. As 
we discuss in more detail later, depends essentially on mi/2, and coannihilation increases 
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by a factor ~ 2 the cosmological upper limit on mi/2 to ~ 1400 GeV, allowing < 600 GeV. 
At this upper limit on mi/2, there is a unique allowed value of mo ~ 350 GeV, but for lower 
values of mi/2 the width of the allowed range of mo expands, reaching 50 < mo ^150 GeV 
when mi/2 ~ 200 GeV. At this value of mi/2 and scanning across the cosmological range in 
mo, we find that m^ ~ 80 GeV, with a small variation by ~ 0.4 GeV. These numbers are 
not very sensitive to tan/5 in the range from 3 to 10 studied in |15| and here, nor are they 
very sensitive to the chosen value of A. 

The lower limit on mi/2 and hence depends on the sparticle search limits provided 
by LEP 1^. The most essential of these for our current purposes are those provided by 
the experimental lower limits on the lighter chargino mass m^± and the lighter scalar Higgs 



mass rriHo 



A lower limit m^± > 95 GeV was assumed in |15]: unsuccessful chargino searches 



during higher-energy runs of LEP have now increased this lower limit to m^± > 100 GeV [21 



which does not reduce very much the range allowed in |15|. 

The impact of the recently- improved lower limits on the Higgs mass [21| is potentially 
more significant, particularly for tan/? = 3, as displayed in Figs. 6 and 7 of |1^. The 
present experimental lower limit for tan/3 = 3 is probably rriHi > 105 to 109 GeV [21|. 
The contours shown in Figs. 6a,b and 7b of were not calculated with the most 
recent two-loop MSSM code so we take the ttih^ = 100 GeV lines in [|I^] as indicative 
constraints. These correspond to mi/2 ~ 340(720) GeV for fi > (<)0, corresponding in turn 
to m^ > 140(310) GeV. On the other hand, for tan/? = 10, the LEP lower limit on 
is considerably weaker than 100 GeV |^], and hence does not constrain significantly the 



allowed parameter space, as seen in Figs. 6c, d and 7c of [|15 



We note in passing that requiring our present electroweak vacuum to be stable against 
transitions to a lower-energy state in which electromagnetic charge and colour are broken 



(CCB) would divide the parameter regions allowed in |T5| into two parts: one at large 
mi/2 and the other at small mi/2 and relatively large mo. We do not implement the CCB 
constraint in our analysis, since it may be considered optional. Nor do we implement any 
constraint due to the observed rate of 6 ^ 57 decay ||2^, but it is well known that this 
reduces very substantially the parameter space allowed for /i < 0. 



5 Results 

As discussed above, we scan the cosmologically preferred set of parameters which yield 
0.1 < Q^h"^ < 0.3 and are consistent with the recent LEP accelerator bounds. For each 
value of tan (3 and sign of /i, we vary mi/2 and mo over all the allowed range. As default, we 
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choose Aq = —nil 12 in most of our computations. Then, using the hadronic inputs described 
in section 3, we compute separately the spin- dependent and scalar contributions from the 
a2 and as coefficients, respectively, to the elastic scattering of LSPs on both protons and 
neutrons. 

In Figure 1, we show the resulting spin- dependent elastic cross section as a function of 
the LSP mass, m^^,. Although it is barely discernible, the thicknesses of the central curves in 
the panels show the ranges in the cross section for fixed that are induced by varying itlq. 
At large where coannihilations are important, the range in the allowed values of mo is 
small and particularly little variation in the cross section is expected. The shaded regions in 
this and the following figures show the effects of the uncertainties in the input values of the 
(H). In Fi gure la, for tan/3 = 3,yU < 0, we see at small the effect of a cancellation 
induced by the difference in signs between A„ and A^ ,,- Cancellations are possible for the 
other values of tan/3 and sign of /x, but not in the preferred range of mi/2 and mo used here. 
Aside from the cancellation, the spin-dependent cross section peaks at about 10~^ pb and 
drops rapidly as increases. 

In Figure 2, we show the corresponding result for the scalar cross section, based on a^. 
As in Figure 1, the thickness of the central curve reflects the range in mo sampled. The 
shaded region now corresponds to the uncertainties in the inputs given in (|r^). The scalar 
cross section is, in general, more sensitive to the sign of ^ than is the spin-dependent cross 
section. Notice that, in Figure 2c for tan/5 = 10 and /i < 0, there is another cancellation. 
In this case, Higgs exchange is dominant in as. We flrst note that, for yU < 0, both Z^s 
and are negative, as is the Higgs mixing angle a. Inserting the deflnitions of 5H(2i), we 
see that there is a potential cancellation of the Higgs contribution to as for both up-type 
and down-type quarks. Whilst there is such a cancellation for the down-type terms, which 
change from positive to negative as one increases m^, such a cancellation does not occur 
for the up-type terms, which remain negative in the region of parameters we consider. The 
cancellation that is apparent in the flgure is due to the cancellation in as between the up-type 
contribution (which is negative) and the down-type contribution, which is initially positive 
but decreasing, eventually becoming negative as we increase m^. 

In Figure 3, we show the effects of varying some of the input assumptions made earlier. 
For example, when the assumed uncertainty in y is taken to be ±0.2, we get a thicker shaded 
region, as shown for tan/5 = 3,/i < in Figure 3a. In Figure 3b, we give one example of 
the cross section for the elastic scattering of neutralinos on neutrons. This particular case 
was chosen because it displays the largest difference between the neutron and proton cross 
sections among those tested. As one can see, our results for the LSP scattering on neutrons 
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Figure 1: The spin- dependent cross section for the elastic scattering of neutralinos on 



protons as a function of the LSP mass. The central curves are based on the inputs ^^), and 
their thicknesses are related to the spreads in the allowed values ofrriQ. The shaded regions 



correspond to the uncertainties in the hadronic inputs (24) 
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Figure 2: The spin-independent scalar cross section for the elastic scattering of neutralinos 
on protons as a function of the LSP mass. The central curves are based on the inputs 1^1^, 
their thicknesses are again related to the spread in the allowed values of m^, and the shaded 
regions now correspond to the uncertainties in the hadronic inputs (ff^. 
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and protons are almost identical. Similarly, the effects of changing Aq are also relatively 
small, as illustrated by two cases with = in Figures 3c and 3d. In the latter example, 
there is almost a factor of 2 difference at higher values of m^, which is due to yet another 
cancellation, this time between the squark-exchange and Z-exchange terms in a2u- 

Finally, we show in Figure 4 compilations of our results for the spin-dependent and 
-independent cross sections, compared with current and projected experimental limits ob- 
tained from [^|. The shaded region in panel (a) is the union of the shaded regions in Figure 



1, and the shaded region in panel (b) is the union of the shaded regions in Figure 2. 

6 Discussion 



As seen in Figure 4, the present experimental upper limit on the spin-independent part of 



the elastic scattering of the LSP on a nucleon is around 10~ pb for 50 GeV < < 100 GeV. 
On the other hand, the maximum scalar cross section we find is around 10~^ pb, which is 
attained for ~ 50 GeV. This means that present experiments searching directly for 
supersymmetric dark matter are far from constraining the parameter space of our baseline 
theoretical framework, in which LEP constraints are applied to MSSM models with universal 
supergravity-inspired soft supersymmetry-breaking parameters mi/2,fno- 

The literature contains predictions for the elastic LSP-nucleus scattering rates that vary 
considerably, with some estimates lying considerably higher than ours P, [l^. There are 



various ways in which such differences might arise, of which we mention a few here. We 
have imposed the requirement that the LSP relic density lie in the favoured range 0.1 < 
^x^^ — 0-3) whereas other calculations often include models with lower relic densities. Such 
models would normally have larger XX annihilation cross sections, and correspondingly larger 
elastic scattering cross sections. Hence the predicted scattering rates would be larger, if the 
conventional halo density p ~ 0.3 GeV/cm^ is assumed for the LSP However, we 

believe this assumption is unreasonable: if not all the total cold dark matter density flcDM 
is composed of LSPs, the density of LSPs in the halo should be reduced by the corresponding 
factor Q^/QcDM- 

Other possible differences may arise in the treatment of the LEP constraints: we find it 
to be almost excluded that the LSP be Higgsino-like 0, even if the assumptions of universal 
soft supersymmetry breaking are relaxed, and Higgsino dark matter is certainly excluded 
if universality is assumed, as is the case here. In addition to the LEP constraints, this is 
because the value of fi is predicted as a function of mi/2 and mo, placing the LSP firmly in 
the Bino-like region. The same considerations exclude an LSP with mixed Higgsino/gaugino 
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Figure 3: As in Figures 1 and 2, hut now illustrating a) the effect of enlarging the uncertainty 
in y to ±0.2, b) the spin- dependent cross section for the elastic scattering of neutralinos on 
neutrons, c) and d) the effects of setting = for < and two choices of taxi 13. 
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Figure 4: Summary plot compiling our predictions in the constrained MSSM for (a) the 
spin- dependent and (b) the spin-independent elastic LSP scattering cross sections on protons, 
compared in each case with current and projected limits on the scattering cross sections, as 



obtained from 



content. 

The prediction of /i may be circumvented by postulating non-universality for the soft 
scalar supersymmetry-breaking parameters in the Higgs sector, which might have appeared 



to resurrect the possibility of a Higgsino-like LSP However, such a possibility goes 

beyond the universality framework adopted here, and, moreover, the LEP constraints now 
appear to exclude this possibility |0, as mentioned above. 

ra, 0, IT8 



use 



There are no differences between the effective Lagrangians we and others 
to describe the four-fermion x~ 1 interaction that determines the elastic x~P-,^ scattering 
cross sections. We have found differences of detail between our and other treatments of 
the hadronic matrix elements of the scalar and axial-current qq operators appearing in this 



Lagrangian [^, but this is not responsible for any big difference between the results. 

We should not want our experimental colleagues to be too downcast by the long road 
they appear to have to cover in order to probe the minimal universal MSSM framework 
utilized here. For example, there are surely some supersymmetric models that predict larger 
scattering rates. However, we think it best to have in mind a plausible and realistic target 
sensitivity, which is what our universal framework and implementation of the LEP and 
cosmological constraints provide. Our results also have the merit of being relatively specific: 
as seen in Figure 4, the elastic scattering cross sections we predict for any given value of the 
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LSP mass lie in a comparatively narrow band. As discussed earlier, this is essentially 
because the LSP is always mainly Bino-like in our framework, so its couplings do not depend 
greatly on other MSSM parameters such as mo- The principal causes of broadening are the 
uncertainties in the hadronic inputs and the possibilities of cancellations that may reduce 
the cross sections for some specific values of the constrained MSSM parameters. 

This tight correlation we find between the LSP mass and its elastic scattering rate means 
that future experiments should be able to phrase their sensitivities directly in terms of 
the LSP mass in the universal supergravity-inspired version of the MSSM. For example, our 



results suggest that the proposed Genius experiment ||28[ would be sensitive to < 100 GeV 
for almost all MSSM parameter choices in Figure 5b. More optimistically, if/when a signal 
is observed, its plausibility would be enhanced if its recoil spectrum was correlated with the 
rate in the manner suggested by Figure 5. Thus our analysis provides experiments with an 
additional tool that may assist in the extraction of a signal that might be significantly smaller 
than they could have hoped. In any case, the importance of the search for supersymmetric 
matter remains unchanged, and there are still several years before the LHC comes into 
operation, so these experiments still have both motivation and opportunity. 
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